arXiv:1506.05334v2 [hep-ph] 23 Jan 2016 


Cosmic axion background propagation in galaxies 


Francesca V. Day a 

°Rudolf Peierls Centre for Theoretical Physics, University of Oxford, 
1 Keble Road, Oxford, 0X1 3NP, United Kingdom 


Abstract 

Many extensions of the Standard Model include axions or axion-like particles 
(ALPs). Here we study ALP to photon conversion in the magnetic field of the 
Milky Way and starburst galaxies. By modelling the effects of the coherent and 
random magnetic fields, the warm ionized medium and the warm neutral medium 
on the conversion process, we simulate maps of the conversion probability across 
the sky for a range of ALP energies. In particular, we consider a diffuse cosmic 
ALP background (CAB) analogous to the CMB, whose existence is suggested by 
string models of inflation. ALP-photon conversion of a CAB in the magnetic fields 
of galaxy clusters has been proposed as an explanation of the cluster soft X-ray 
excess. We therefore study the phenomenology and expected photon signal of CAB 
propagation in the Milky Way. We find that, for the CAB parameters required to 
explain the cluster soft X-ray excess, the photon flux from ALP-photon conversion 
in the Milky Way would be unobservably small. The ALP-photon conversion prob¬ 
ability in galaxy clusters is 3 orders of magnitude higher than that in the Milky 
Way. Furthermore, the morphology of the unresolved cosmic X-ray background is 
incompatible with a significant component from ALP-photon conversion. We also 
consider ALP-photon conversion in starburst galaxies, which host much higher 
magnetic fields. By considering the clumpy structure of the galactic plasma, we 
find that conversion probabilities comparable to those in clusters may be possible 
in starburst galaxies. 
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1. Introduction 

Axions and axion-like particles (ALPs) arise in many extensions of the Stan¬ 
dard Model as pseudo-Nambu-Goldstone bosons of broken symmetries. A generic 
ALP is an ultra-light pseudo-scalar singlet under the Standard Model gauge group. 
Throughout this work we will consider massless ALPs with no coupling to QCD. In 
the low energy effective field theory, an explicit mass term in the Lagrangian is for¬ 
bidden by a shift symmetry a{x) a{x) -(-constant. We expect non-renormalizable 
couplings between the ALP and the Standard Model suppressed by the high scale 
M. In this work we explore the phenomenology of the dimension five 077 coupling. 
In addition to the Standard Model Lagrangian, we therefore have: 

£alp = ^ad^a + ~^aFF = ^ad^a + -^E • B. (1) 
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The term C D jjE • B leads to ALP-photon interconversion in the presence of 
a background magnetic field. The mass scale M is model dependent and so is a 
priori undetermined. Empirical limits on M may be derived from astrophysical 
observations and from axion search experiments, as reviewed in |T]. For low mass 
ALPs ( m a < 10~ 10 eV), the strongest bounds on M arise from observations of the 
SN1987a supernova in the Large Magellanic Cloud. In ALP extensions of the Stan¬ 
dard Model, we would expect an ALP burst coincident with the neutrino burst. 
This ALP burst would be observable as a gamma ray flux following ALP-photon 
conversion in the Milky Way magnetic field. The non-observation of such a gamma 
ray flux leads to the bound M > 2 x 10 11 GeV [2HIj. 

A primordially generated, thermally produced cosmic ALP background (CAB), 
analogous to the CMB, is a natural prediction of string theory models of inflation 
[5], The CAB has a quasi-thermal energy spectrum that is red shifted to soft X- 
ray energies today. The constituent ALPs act as dark radiation - extra relativistic 
degrees of freedom conventionally parametrised by the equivalent excess in the 
number of neutrino species AN e Q. Current measurements of AN e g are consistent 
both with zero and with a significant dark radiation component [6]. The ALP 
number density in the CAB between energies E and E + dE is: 


dN (E) = AX (E) dE, (2) 

where X (E) is the shape of the CAB energy spectrum and A its normalisation. 
The spectral shape is predicted by the general string inflation scenario described 
in [5], and may be found by numerically solving the Friedmann equations for ALP 
production and redshift. The resulting spectrum is fit well by the function 


X(E) = E q e~ aEr . (3) 

The constants q , a and r are found by fitting equation ([3]) to a numerical solution 
of the equations of motion, and in general depend on the mean ALP energy Ecab- 
In a typical string inflation model, Ecab ~ 0(100 eV). The overall normalisation 
of the spectrum is model dependent but may be measured by its contribution to 
AiVeff. We will therefore find the normalisation constant A by setting the CAB 
contribution to AN e g. This is related to the CAB energy density by: 


PCAB = A JVgff- 

O 


(n) ’ 


The flux d<3? a (E) of ALPs with energies between E and E + dE is then: 


(4) 


d&a (E) = dN (E) (5) 

so, 

— =AX(E) C . (6) 

dE K ’4 V 7 

The predicted spectrum of the CAB background for Ecab = 200 eV and AA^ = 

0.5 is shown in figure [lj In this case, the parameters in equation ([3]) are found to 

be q = 0.62, r = 1.5, o = 2.6 x lO^eV” 1 ' 5 . 
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For sufficiently large jj and CAB flux, ALP-photon conversion offers the pos¬ 
sibility of detecting a CAB as an excess of soft X-ray photons from environments 
with a sufficiently strong and coherent magnetic field [7]. A natural place to search 
for this effect is in galaxy clusters, which host 1 — 10 /xG fields over Mpc distances. 
Furthermore, there is a long standing excess in the soft X-ray (E < 400 eV) flux 
observed from galaxy clusters, above the predicted thermal emission from the intra¬ 
cluster medium. It was suggested in [7j that CAB to photon conversion in galaxy 
clusters could be the source of this soft X-ray excess. Detailed simulations of this 
process have been carried out for the Coma [HI [9j, A665, A2199 and A2255 [IU] 
galaxy clusters. These show that CAB to photon conversion can consistently ex¬ 
plain the observed excess in Coma, A2199 and A2255 as well as the non-observation 
of an excess in A665, within astrophysical uncertainties. 

While by no means conclusive, this hint of new physics motivates studying 
the consequences of a CAB in other astrophysical systems. Galaxies also host 
magnetic fields and are therefore potential ALP to photon converters, as discussed 
in mmmm- Note that we come to qualitatively different conclusions than those 
in JT5]. This is discussed further in the Appendix. ALPs from a CAB may convert 
to X-ray photons in the Milky Way. This would contribute to the unresolved 
cosmic X-ray background - the diffuse X-ray intensity observed across the sky 
after subtracting the integrated emission from all detected point sources. Within 
standard physics, the unresolved cosmic X-ray background could arise from the 
Local Bubble and the warm-hot intergalactic medium mm • There is also room 
for more exotic contributions, such as decaying dark matter or the CAB considered 
here. However, unlike in the case of the galaxy cluster soft X-ray excess, there is no 
problem explaining the cosmic X-ray background within the framework of standard 
physics. The possibility that the cosmic X-ray background is related to conversion 
of a CAB to photons in the Milky Way’s magnetic field was first considered in 

mm- 

The magnetic field in a starburst galaxy (a galaxy with a very high rate of star 
formation) is typically an order of magnitude higher than that in the Milky Way, 
suggesting a substantially higher rate of ALP-photon conversion. We therefore also 
estimate the ALP to photon conversion probability in starburst galaxies. In both 
cases, the warm ionized and neutral gas in the galaxy also plays a significant role 
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in determining the ALP to photon conversion probability, as discussed in sections 
2 and 5. 

In this paper, we will discuss the phenomenology and potential observational 
consequences of a CAB’s passage through the Milky Way. We will also consider 
ALP to photon conversion in the high magnetic field, high plasma density environ¬ 
ment of starburst galaxies. In section 2 we consider in more detail the propagation 
of the ALP-photon vector in galactic environments, focusing on the effect of the 
magnetic field, the warm ionised medium and the warm neutral medium. In section 
3, we describe our model of the Milky Way environment. In section 4, we present 
and discuss our results for ALP-photon conversion in the Milky Way. In section 
5, we discuss some caveats and additional relevant effects. In particular we derive 
the conditions under which the dumpiness of the warm ionized gas in galaxies 
becomes relevant for ALP-photon conversion. In section 6 we apply this to esti¬ 
mate the ALP-photon conversion probability in starburst galaxies. We conclude in 
section 7. 


2. ALP-photon conversion 


The ALP-photon coupling is suppressed by an energy scale M much larger than 
the physical energies involved. It is therefore sufficient to simulate ALP-photon 
conversion using the classical equation of motion derived from (JT]) , and neglecting 
higher dimension terms. We further assume that the ALP wavelength is much 
shorter than the scale over which its environment changes, allowing us to linearise 
the equations of motion. This condition is abundantly satisfied for X-ray energy 
ALPs in astrophysical environments. The ALP-photon equations of motion are 
conveniently described by combining the ALP with the two photon polarizations 
in an ALP-photon vector. The linearised equation of motion for an ALP-photon 
vector of energy cj is then: 


l 

( A 7 
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ALP-photon mixing 
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^7 ai 


— id 
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= 0 . 


(7) 


rections perpendicular to the direction of travel. Ap describes Faraday rotation 
between the two photon polarizations. This effect is not relevant to ALP-photon 
conversion, and so we set = 0. The photon components are given an effective 
mass by their interactions with free electrons in the surrounding medium. This 
effective photon mass is equal to the plasma frequency - the frequency of charge 


density oscillations in the surrounding plasma. This is given by u p i = ( 47 tq: 
where n e is the free electron density. We then have : 


A _ ~ u li _ 4:iran e 
2u ; 2 u>m e 

As we do not measure the photon polarization, we simply add the conversion prob¬ 
abilities for each polarization. For an initially pure ALP state, in our semi-classical 
approximation the conversion probability after a distance L is: 
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Pa^(L) = I (1,0,0|/(L)> I 2 + I (0,1, 0|/(L)> I 2 , (9) 

where | f(L)) is the final state after a distance L as determined by equation (ITJ). 

b 2 +b 2 

The ALP to photon conversion probability P a -> 7 is proportional to X M $ y in the 
limit jj <C 1. A non-zero electron density in the propagation environment gives 
an effective mass to the photon, causing decoherence between the ALP and photon 
components and hence suppressing P a _> 7 . 

For constant electron density and magnetic field, there is an analytic solution for 
the conversion probability. We identify two angles associated with the propagation: 


tan ( 2 0) = 2.8 x 10 -3 x 

A = 0.053 x 


n, 


10 3 cm 3 

1 keV 


n e 


) (^) (t£v) ( 


w \ / 10 i3 GeV\ 

M ) 


10 3 cm 3 


UJ 


L 


1 kpc 


For a single domain of length L, the conversion probability is then 


P(a —> 7 ) = sin 2 (29) sin 2 


A 


cos 29 


,( 10 ) 

( 11 ) 


( 12 ) 


In a more general case, for P a _>. 7 < 1, we find: 


Pa^(L)=y\ [ L dze^A^aiiz) 

• Jo 


I=x,y 


(13) 


where, 

<P( Z ) = f dz'A^z') = --jj J dz'ujpi(z'). (14) 

As shown in equation ([ 8 ]), A y (z) oc n e , and so the electron density has the effect of 
rotating the probability amplitudes ( 1 , 0 , 0 |/(L)) and ( 0 , 1 , 0 |/(L)) in the complex 
plane as L increases, suppressing the efficacy of the magnetic field in increasing the 
conversion probability over increasing distances. 


The Milky Way is almost opaque to low energy X-rays due to photoelectric 
absorption from the warm neutral medium. We capture this effect in our equation 
of motion for the ALP-photon vector by including a damping parameter r(z) that 
describes the attenuation of the photon components. The new equation of motion 
no longer describes a closed quantum system - the Hamiltonian for the ALP-photon 
vector alone is no longer Hermitian. We therefore use a density matrix formalism: 

/ A 7 0 A 7ax \ ( i\ 0 0 \ 

H= I 0 A 7 A iay - 0 7 ^ 0 , (15) 

V A 7ax A iay 0 / V 0 0 0 / 

( I 7.) \ 

P = I 72/) ® ( | 7z> IT y) I a) ) , (16) 

V I a) J 

p(z) = e~ iHz p(0)e iHU . (17) 
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To simulate the conversion probabilities, we discretize each line of sight into do¬ 
mains of length 6z: 


Pk 


„—iHu5z „ iH}5z 
e Pk-ie k , 


(18) 


where p^ is the density matrix in the fcth domain and is the Hamiltonian defined 
using the magnetic field, electron density and neutral hydrogen density in the centre 
of the kth domain. 


3. The Milky Way environment 

Three properties of the Milky Way’s interstellar medium are relevant to ALP- 
photon conversion - the magnetic field, the free electron density provided by the 
warm ionized medium and the opacity to X-rays provided by the warm neutral 
medium. In this section, we describe our model for each of these components. We 
leave a discussion of various caveats to and justifications of this model, in particular 
the dumpiness of the electron density, to section 5. 


3.1. Magnetic field 

We use the recent model by Jansson and Farrar [201121] . based on 40,000 extra- 
galactic Faraday rotation measures. The magnetic field is the sum of three compo¬ 
nents - the coherent field, the random field and the striated held. The coherent held 
has large scale structure on the scale of the Milky Way with typical held strengths 
of a few pG. The coherent held is modelled as the sum of a disc held, which follows 
the spiral arms of the Milky Way; a halo held above and below the disc; and an 
‘X held’ which points out of the plane of the Milky Way. The radial extent of 
the halo held is much greater in the South of the galaxy than in the North. The 
coherent held model in [20] artificially excludes the central 1 kpc of the Milky Way. 
We therefore augment the model with a 5 p.G radially constant poloidal held with 
vertical scale height 1 kpc in the central 1 kpc only. A full sky map of the average 
coherent held is shown in figure [2] 

The random held has a set magnitude with a disc and halo component, but its 
direction is randomized with a coherence length of 0(100pc), the typical size of 
a supernova outhow. The magnitude of the random held is typically a few times 
higher than that of the coherent held. The striated held has a magnitude 1.2 
times that of the coherent held with its sign randomized on coherence scales of 
0(100pc). We see from equation (13) that the conversion probability increases 
with the coherence length of the magnetic held. Indeed, for the majority of the 
Milky Way the coherent held gives the dominant contribution to ALP-photon con¬ 
version. The exception to this is in the disc of the Milky Way, where the random 
held is 0 (10 pG) whereas the coherent held is 0 (1 pG). Additionally, the coherent 
held often reverses sign between the spiral arms, decreasing its coherence length 
in the disc. We therefore use all three held components in modelling ALP-photon 
conversion in the Milky Way. The random and striated helds are implemented with 
respect to each line of sight - the direction and sign respectively are randomised 
every 100 pc along each ALP-photon path separately. This simple implementation 
clearly does not give a realistic picture of the random and striated helds across the 
Milky Way, but is adequate for modelling their effects on ALP-photon conversion. 


6 



Average Transverse Field/pG 



Figure 2: The average coherent transverse magnetic field for lines of sight starting 20 kpc from the 
Earth. Galactic latitude increases vertically in the plot, while galactic longitude increases to the 
right. The centre of the plot corresponds to a line of sight in the direction of the galactic centre. 


3.2. Electron density 

As described above, the photon gains an effective mass through interactions 
with surrounding free electrons. This mass suppresses ALP-photon conversion, as 
shown in equations ([8]) and (13). We use the thin and thick disc components of the 
NE2001 [22] model of the Milky Way electron density: 


f/t.hick (f) 


COS (34^) 

COS (Mkfc) 


H (17kpc 


5thin (r) = e 


r —3.7 kpc N ^ 
1.8 kpc J 


r ), 


n t hick (r, z) = 0.035 cm 3 fif th ick (r) sech 2 

^thin (r, z ) = 0.09 cm” 3 c/ thin (r) sech 2 ^ 
n e (r, z) = n t hick (r, z) + n thin (r, z ), 


0.95 kpc J 

z 


0.14 kpc J ’ 


(19) 


where (r, z) are cylindrical polar coordinates centred at the galactic centre, Rq = 
8.5 kpc is the distance to the Sun and H (x) is the Heaviside step function. This 
model predicts unphysically low electron densities at large radii. While this is 
not important for many astrophysical phenomena, which depend only on line of 
sight integrals of n e , it can have a large effect on P a ^. 7 . We therefore enforce a 
minimum electron density of n m i n = 10 -7 cm -3 , approximately the electron density 
of inter-galactic space. 


3.3. Photoelectric absorption by the warm neutral medium 

As explained in section 2, we model photoelectric absorption with the damping 
parameter T, which describes the attenuation of the photon component of the ALP- 
photon vector. This is conventionally parameterzied by the effective cross section 
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with respect to neutral hydrogen, so that T (w) = <r e ff (w) (nm + 2n# 2 ), where 
nm + 2 njj 2 is the density of neutral hydrogen. (HI refers to atomic hydrogen 
and H 2 to molecular hydrogen.) Photoelectric absorption by heavier elements 
(which is dominant for ui > 1 keV) is included in the effective cross section <r e fr (w) 
by assuming solar abundances for the relative densities of hydrogen and heavier 
elements. We use effective cross section values from [23] - we note in particular that 
cr e fj (w) is highly energy dependent, ranging from cr e fr(100eV) = 5.7 X 10 _2 ° cm 2 
to <7 e ff (2keV) = 4.3 x 10 _23 cm 2 . We use the neutral hydrogen densities given in 

m ■■ 


nm 


0.32 cm 3 exp ^ 18 . 2 4k P c o.52k P c) 
0 , 


if r > 2.75 kpc 
otherwise 


nH 2 = 4.06 cm 3 exp I — 


1 z 1 


2.57 kpc 0.08 kpc 


( 20 ) 

( 21 ) 


4. The Milky Way: Results and Discussion 


We now apply equations (15) and (18) to the ALP-photon conversion in the 
Milky Way, as described above. 


4-1- Conversion Probabilities 

The ALP to photon conversion probability across the Milky Way is relevant to 
a range of effects, probing different ALP energies. Our focus in this paper is on 
the propagation of a cosmic ALP background at soft X-ray energies converting to 
photons in the Milky Way. The Milky Way ALP to photon conversion probability 
at a range of soft X-ray energies is shown in figures and [5} We see that 
at w = 200 eV, P a _> 7 is heavily suppressed by photoelectric absorption but this 
suppression is not significant at u = 500 eV, where P a -ty inherits the geometry of 
the galactic magnetic field. 


Log 10 (Conversion Probability) 



Figure 3: The ALP to photon conversion probability in the Milky Way for ALP energy u> = 200 eV 
and M = 10 13 GeV. 
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Log 10 (Conversion Probability) 



Figure 4: The ALP to photon conversion probability in the Milky Way for ALP energy uj = 500 eV 
and M = 10 13 GeV. 



Figure 5: The ALP to photon conversion probability in the Milky Way for ALP energy ui = 1500 eV 
and M = 10 13 GeV. 


It has been suggested that the 3.5 keV line recently observed in galaxy clusters 
and the Andromeda galaxy may arise from dark matter decay to ALPs followed 
by ALP-photon conversion in astrophysical magnetic fields PU OH [251 [26]. This 
scenario fits the morphology of the 3.5 keV line flux in galaxy clusters and predicted 
its non-observation in a stacked sample of galaxies. In this case, the conversion 
probability for 3.5 keV ALPs in the Milky Way, shown in figure [6] is required. As 
shown in jl'6j . the conversion probability is too low to expect an observable 3.5 keV 
line flux in the Milky Way halo. 

ALP-photon conversion has also been suggested as an explanation for the 
anomalous transparency of the universe to gamma rays EMU EMoj- In this 
scenario, gamma rays emitted by distant blazars convert to ALPs in the magnetic 
field of the host galaxy or in the intergalactic magnetic held, and then reconvert to 
photons in the intergalactic or Milky Way magnetic held. In this way, gamma ray 
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Log 10 (Conversion Probability) 



Figure 6: The ALP to photon conversion probability in the Milky Way for ALP energy uj = 3.5 keV 
and M = 10 13 GeV. 


photons are able to avoid scattering from electrons in intergalactic space. In this 
case, the conversion probability for gamma ray energy ALPs, shown in figure [7j is 
key. This conversion probability is also used in calculating the bounds on M from 
SN1987a. Comparing figures ®® and [7| we see that the conversion probability 
in the Milky Way saturates at P a _>. 7 ~ 10~ 6 for M = 10 13 GeV. This behaviour 
can be seen in the single domain formula (equation (| 12|) ), where for u —> oo, the 
analytic formula for reduces to: 


lim P a _= sin 2 

UJ—>00 



( 22 ) 


Log, 0 (Conversion Probability) 



I 



Figure 7: The ALP to photon conversion probability in the Milky Way for ALP energy u> = 1 TeV 
and M = 10 13 GeV. 
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4-2. Application to a cosmic ALP background 


A cosmic ALP background would convert to photons in astrophysical magnetic 
fields leading to a potentially observable soft X-ray flux. This effect has been sug¬ 
gested as the source of the soft X-ray excess in galaxy clusters. The ALP to photon 
conversion probability in the Milky Way is around 3 orders of magnitude lower than 
that in galaxy clusters, primarily due to the Milky Way’s smaller size. We therefore 
do not expect such a strong signal from CAB to photon conversion in the Milky 
Way. Any extra soft X-ray photons generated from a CAB’s passage through the 
Milky Way would contribute to the unresolved cosmic X-ray background - the dif¬ 
fuse soft X-ray flux remaining after subtracting the flux from known point sources. 

The photon flux from a CAB in a solid angle 12 is given by: 


d®n 

dE 


1 

2n 



dtl'P (n',E) AX (E) 


c 

4' 


(23) 


For example, for a central energy Eqab = 200 eV, the predicted photon fluxes for 
disc and halo pointings are shown in figures [8] and [9j We normalise to AA r e ff = 0.5 
to allow easy comparison with the galaxy cluster fluxes simulated in |8j. The pre¬ 
dicted CAB signal scales linearly with the CAB contribution to AN e g. We notice 
that the shape of the spectrum is dramatically altered from the CAB spectrum 
shown in figure [TJ as the conversion probability at low energies is dramatically 
suppressed by photo-electric absorption. The spectrum is further influenced by 
the energy dependence of the conversion probability even in the absence of ab¬ 
sorption. For example, for both pointings we see oscillations in the predicted 
flux on top of the overall shape of the spectrum. These can be understood by 
considering the analytic solution in equation (12), which approximates the qual¬ 
itative features of the solution in the non-homogeneous case simulated here. In 
particular, we expect to see local maxima in the conversion probability whenever 

A = 0.053 x ( lo-a'cm-s ) (l¥^) = for odd integer N. These correspond 

to the oscillations seen in figures [8] and [9] and are a distinctive feature of a photon 
flux arising from ALP to photon conversion in a sufficiently high electron density 
environment (so that A > 1). The flux from the Milky Way centre (figure is 
lower and peaks at higher energies that that from due Galactic North (figure (9| due 
to the higher warm neutral medium column density towards the Galactic centre, 
leading to a greater effect from photo-electric absorption. Note that the detailed 
shape of the Milky Way centre spectrum is highly dependent on the realisation of 
the strong random magnetic field in the disc. 
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Photon Flux/eV 1 s Vm 2 sr 1 



Figure 8: The predicted photon flux from the Milky Way centre from CAB to photon conversion 
with Ecab = 200 eV and M = 10 13 GeV 


Photon Flux/eV 1 s 1 cm 2 sr 1 



Figure 9: The predicted photon flux from due Galactic North from CAB to photon conversion 
with Ecab = 200 eV and M = 10 13 GeV 


Full sky maps of the cosmic X-ray background were observed by the ROSAT 
satellite m- We now calculate the predicted flux from a CAB in the ROSAT 1/4 
keV and 3/4 keV bands. We use equation (23) with the conversion probabilities 
calculated using equation (18). As shown in iSHTOj, natural CAB parameters to 
explain the cluster soft excess are Ecab = 200 eV and M = 5 x 10 12 GeV for 
AlVeff = 0.5. We plot the predicted ROSAT signals for these parameters as full sky 
maps in figures [To] and 11 Comparing with [31], we find that the predicted CAB 
signal is over 3 orders of magnitude smaller than the signal observed by ROSAT. 
The soft X-ray excess can also be explained with a lower Eqab and lower M - in 
this case the signal in the Milky Way is even lower due to the higher photo-electric 
absorption at lower energies. 
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Log 10 [Expected Flux/cm 2 s ^ 



Figure 10: The predicted photon flux in the ROSAT 1/4 keV band for Ecab = 200 eV, M = 
5 x 10 12 GeV and AN e s = 0.5 



Figure 11: The predicted photon flux in the ROSAT 3/4 keV band for Ecab = 200 eV, M = 
5 x 10 12 GeV and AN e g = 0.5 

We might wonder if a CAB with different parameters could contribute signifi¬ 
cantly to the cosmic X-ray background (quite possibly by ignoring the problem of 
overproduction in clusters). The cosmic X-ray background is most clearly seen in 
the Chandra Deep Field (CDF) observations |T8] . The observed fluxes, ALP to pho¬ 
ton conversion probabilities and predicted fluxes from CAB to photon conversion 
(with the parameters used above) for the CDF-South and CDF-North observations 
are shown in table 1. To simulate the conversion probabilities here, we did not in¬ 
clude the effects of photo-electric absorption, as the CDF pointings are chosen for 
their low warm neutral medium column density. For any CAB parameter values we 
expect the CAB signal from CDF-North to be O (10) times lower than that from 
CDF-South. However, the cosmic X-ray background intensities from these observa¬ 
tions are the same within their errors. Therefore the possibility of a CAB forming 
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the dominant part of the cosmic X-ray background is excluded by the North-South 
asymmetry of the Milky Way magnetic field. 



(M) 

P 2—>7 

Predicted Flux 

Observed Flux 

CDF-South 

(- 22 . 8 °, 161°) 

2.4 x 1(T 6 

2.4 x 10 ~ 15 

( 1.1 ± 0 . 2 ) x 10~ 12 

CDF-North 

(54.8°, 125.9°) 

2.5 x 1(T 7 

2.0 x 10~ 16 

(9 ±3) x 10 ” 13 


Table 1: The ALP to photon conversion probabilities P a -> 7 (averaged over the energy range), 
predicted CAB fluxes and observed cosmic X-ray background fluxes after point source subtraction 
for the Chandra Deep Field pointings. Fluxes are given in units of ergs cm -2 s _1 deg -2 . We 
also show the galactic latitude b and longitude l of the observations. We use Scab = 200 eV, 
M = 5 x 10 12 GeV and AN e g = 0.5. 


We see that a CAB responsible for the cluster soft X-ray excess would currently 
be unobservable in the Milky Way, and that a CAB cannot contribute significantly 
to the observed unresolved cosmic X-ray background without giving it a North- 
South asymmetry ruled out by observations. Ubiquitous features of a CAB Milky 
Way signal are a prominent North-South asymmetry (as shown in figures 10 and 
[IT] ), and complex spectral shapes from a convolution of the quasi-thermal CAB 
spectrum and the energy dependent conversion probability as shown in figures [ 8 ] 
and [9] In particular, the conversion probability and therefore the predicted flux 
oscillates as the energy increases. 


5. Additional effects 

5.1. Milky Way magnetic field 

We have used a simplistic model for the random and striated fields with a 
single coherence length of 100 pc. In reality, we expect these fields to exhibit a 
range of coherence scales. However, changing the coherence length by a factor of 
10 in either direction only results in a factor of < 2 difference in the full conversion 
probability. Furthermore, we have not considered the field at the very centre of 
the Milky Way. ALP to photon conversion in the Milky Way centre is discussed in 
m in the context of the 3.5 keV line. Here we simply note that estimates of the 
magnetic held in the Galactic centre are highly uncertain, ranging from 10 /iG to 
1 mG. At the upper end of this held range, we could see conversion probabilities in 
the Milky Way centre comparable to those in galaxy clusters, and therefore might 
expect an observable soft X-ray flux from a CAB. However, the high density of the 
warm neutral medium in the Galactic centre would significantly suppress the signal 
at low energies. Furthermore, the galactic centre is a highly complex environment 
so it would be very difficult to pick out a small excess soft X-ray hux. 


5.2. Clumpiness of the interstellar medium 

We recall that the electron density of the surrounding medium suppresses ALP- 
photon conversion by giving an effective mass to the photon, as shown in equations 
0 , ( |13[ ) and (14). The electron density model used to simulate ALP-photon 

conversion describes the smooth, volume averaged electron density. In reality, the 
warm ionized medium (WIM) in galaxies exists in high density clouds with a rather 
low intercloud electron density |32j . This structure is characterized by the filling 
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factor /, the fraction of a line of sight occupied by WIM clouds. In principle, 
by using the clumpy local electron density n e we might predict a different P a _>. 7 
than we would have by naively implementing the smooth volume averaged elec¬ 
tron density n e . To examine the effect of the local electron density distribution, 
we consider the role of the electron density in rotating the probability amplitude 
A(L) = (1,0,0|/(L)) in the complex plane. For simplicity, we consider the case of 
a constant magnetic field in the x direction, so that only x polarized photons are 
produced. The relevant equations are then (see equations (13) and (fl4|)): 



Pa-ry(L) — 

rL 

/ dze^ z) 

Jo 

B x (z) 

2 

(24) 


2 M 

1 

where, 

<P( Z ) = J dz'Aj(z') = --jj 

[ dz’ 

Jo 

2 t i\ 

Upllz ) , 

(25) 

with 

= 47ra—. 

pL m 



(26) 


We see that the angle of turn in the complex plane is given by ip(z), which is linear 
in n e (z). We first note that whether this turning happens continuously or in steps 
does not significantly effect P a ~> 7 . This is demonstrated in figure [12} where we plot 
in the complex plane the probability amplitude A(L) = dze l ' p ( z '> A ia i(z) for a 
propagation distance L = 0 — 1 kpc increasing along the line. We use B± = 1 juG, 
n e = 0.05 cm -3 , cloud size d c = 10 pc, / = 0.1 and u = 500 eV. In the left 
hand plot we use the volume averaged electron density, and in the right hand plot 
implement the electron density in evenly spaced clouds, with an intercloud electron 
density of 10 - ' cm -3 . We see that the overall shape of A(L) is the same in each 
case, although in the volume averaged case the turn is continuous, whereas in the 
right hand plot we see corners (with a very high rate of turn) where there is a 
cloud. As expected, the conversion probability at L = 1 kpc is practically the same 
in each case (4.9 x 10 -10 and 5.2 x 10 -10 respectively for M = 10 13 GeV). 

However, e i<p ^ is also periodic in n e (2), and it is in this periodicity that we 
see the effect of the dumpiness of the WIM. For high electron densities, low filling 
factors and/or low ALP energies it may be that within a single cloud <p(z ) changes 
by > 27r. A cycle of Arg (e i(fi ( z P within a cloud that returns to its starting point does 
not significantly decrease the overall conversion probability. In the regime where 
this phenomenon occurs, the overall large scale turning of A in the complex plane 
is decreased by the organisation of the WIM into clouds. The significant (> 2ir) 
turning within a cloud essentially gives us a ‘free lunch’ - the volume averaged 
electron density is increased, but there is no contribution to the net large scale 
turning, and so the overall conversion probability is not significantly decreased. 
The predicted conversion probability is therefore significantly increased by taking 
into account the cloud structure of the WIM. This effect is demonstrated by figure 
[ 13 } Here we use the same parameters as in figure [l2| but with uj = 100 eV (so that 
A 7 is increased by a factor of 5). We now see the turns within clouds in the right 
hand plot, allowing A(L) to reach larger radii in the complex plane, and increasing 
P a _>. 7 (L) from 3.0 x 10 -11 with a smooth volume averaged WIM to 9.2 x 10 -10 
with a more realistic WIM profile. In spite of the very high electron density within 
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clouds n c = y 1 , the low electron density intercloud regions allow |^4(L)| to grow in 
this regime. 



Figure 12: A(L ) = dze tlp ^ z> A 7a i(z) for L = 0 — 1 kpc increasing along the line. In the left 

hand plot we use the volume averaged electron density, and in the right hand plot implement the 
electron density in evenly spaced clouds, with an intercloud electron density of 10 -7 cm -3 . The 
clouds correspond to the ‘corners’ in the plot. We see that P a -> 7 = \A\ 2 is not significantly effected 
by the presence of clouds. We use B± = 1 [iG, volume average electron density n e = 0.05cm -3 , 
cloud size d c = 10 pc, filling factor / = 0.1, M = 10 13 GeV and co = 500 eV. For the volume 
averaged case, P a -> 7 ( L = lkpc) = 4.9 x 10 -1 °. With clouds, P a -> 7 (L = lkpc) = 5.2 x 10 -10 . 



Figure 13: A(L) = f g L dze^ (z) A lai {z) for L = 0—1 kpc increasing along the line. In the left 
hand plot we use the volume averaged electron density, and in the right hand plot implement the 
electron density in evenly spaced clouds, with an intercloud electron density of 10 -7 cm -3 . The 
clouds correspond to the loops in the plot. We see that P a -> 7 = \A\ 2 is significantly higher when 
the cloud structure is taken into account. We use B± = 1 ^tG, volume average electron density 
n e = 0.05cm -3 , cloud size d c = lOpc, filling factor / = 0.1, M = 10 13 GeV and = lOOeV. 
For the volume averaged case, P a -> 7 (L = lkpc) = 3.0 x 10 -11 . With clouds, P a ^ry (L = lkpc) = 
9.2 x 10 -lt) . 
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We therefore see that the condition for the cloud structure of the WIM to be 
significant is: 


5 = 1.1 x 10 -2 


n. 


10 -2 cm -3 


keV 


u 


( cloud size 
\ 10 pc 


> 2vr 


(27) 


In figure 12 <5 = 1.1, whereas in figure [I3| 6 = 5.5. This condition is almost never 
satisfied in the Milky Way, so in this work we simply use the volume averaged elec¬ 
tron density given in [22] . However, this effect is significant in high electron density 
environments such as starburst galaxies. Furthermore, an analogous effect will op¬ 
erate whenever oscillations are suppressed by an effective mass from astrophysical 
plasmas. 

The warm neutral medium responsible for photoelectric absorption also has 
a cloud-like structure. We find that using a clumpy rather than homogeneous 
warm neutral medium for Milky Way densities only has a significant effect (after 
averaging over cloud positions) on P a _» 7 for u < 200 eV. At these energies, we 
found that photoelectric absorption reduces the expected signal to negligible levels 
in either case. We therefore simply use the volume averaged warm neutral medium 
densities. 


6. Starburst galaxies 

We now consider ALP to photon conversion in starburst galaxies. Starburst 
galaxies host strong magnetic fields of up to O (100 /rG) in the core regions with 
somewhat lower fields in the halo [531 !3lj. making them potentially very good 
ALP to photon converters. However, the fields in starbust galaxies are largely 
turbulent with little or no coherent field. Furthermore, the electron density is 
correspondingly higher at O (100 — 1000 cm' 3 ) m- Naively, we might expect this 
high electron density to substantially suppress P a ^ 7 . However, in this regime 
the cloud structure of the WIM becomes highly significant, as shown in section 
5. The intercloud electron density it also crucial here. For example, consider a 
simplified case with a random field B = 150 fj ,G over a distance of 3 kpc and 
coherent over 100 pc, implemented as described in section 3. We use a volume 
averaged electron density n e = 1000 cm -3 . Using u = 1 keV, M = 5 x 10 12 GeV and 
the constant, volume averaged electron density gives P a ^ 7 ~ 10 -11 , averaging over 
100 instances of the random field. Now let us assume the WIM exists in 1 pc clouds 
(for example as in the starburst galaxy M82 [35]) with filling factor / = 0.1. In 


this case, applying equation (27) we obtain 5 = 1100 and so the presence of clouds 
is of great importance. If we assume an intercloud electron density of 0.1cm -3 
we obtain P a -t 7 ~ 10 -4 - a conversion probability comparable to that in galaxy 
clusters. However, if we assume that the intercloud electron density is 10 cm -3 we 
obtain P a _> 7 ~ 10 -7 , rendering any signal unobservable (although still 4 orders 
of magnitude higher than for a volume averaged electron density). It is therefore 
possible, but certainly not guaranteed, that we might see signals from a CAB in 
some starburst galaxies. Furthermore, at ALP energies E < 500 eV, any signal 
would by highly suppressed by photoelectric absorption. Starbursts might also be 
good observation targets for the 3.5 keV line arising from dark matter decay to 
ALPs discussed in HSUm [251126]. 
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7. Conclusions 


We have simulated ALP to photon conversion probabilities for axion-like parti¬ 
cles propagating through the Milky Way to Earth. We find that the cosmic axion 
background motivated by string models of the early universe and by the cluster soft 
X-ray excess would be entirely unobservable following ALP-photon conversion in 
the Milky Way’s magnetic field. This is due to low conversion probabilities in the 
Milky Way relative to galaxy clusters, as well as the high photoelectric absorption 
cross section for the central CAB energies. Furthermore, the North-South asymme¬ 
try in this magnetic field is not reflected in observations of the unresolved cosmic 
X-ray background, ruling out a significant ALP contribution to the cosmic X-ray 
background. The smaller size of galaxies compared with galaxy clusters make them 
in general poorer targets for observation of ALP-photon conversion. One exception 
might be starburst galaxies, which feature very high magnetic fields and electron 
densities. 

The galactic electron density suppresses conversion by giving an effective mass 
to the photon component, but in such high density environments the detailed lo¬ 
cal structure of the plasma must be considered. We have derived a condition for 
when the cloud structure of a galaxy’s electron density is relevant for ALP-photon 
conversion. We find that in the Milky Way, and other typical spiral and elliptical 
galaxies, the cloud structure is not relevant. However, the cloud structure is rel¬ 
evant in high electron density environments such as starburst galaxies. We found 
that when the cloud structure of the electron density is taken into account, the 
predicted ALP-photon conversion probability in starburst galaxies is increased by 
up to 8 orders of magnitude, depending on the assumed intercloud electron density. 
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Appendix A. Comparison with previous results 

Our main conclusion - that a cosmic axion background cannot contribute sig¬ 
nificantly to the unresolved cosmic X-ray background - differs to that found by 
Fairbairn in [15] (hereafter F2013). Here we explain further why our results differ. 
The key elements are: 

• Consistent conversion of electromagnetic quantities from natural units to SI 
units 

• Modeling of photoelectric absorption using the density matrix formalism 

• The ALP spectrum 
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• The morphology of the unresolved X-ray background 
We now consider each of these factors in turn. 


Appendix A.l. Conversion of B from natural to SI units 

To compute conversion probabilities from equation [7] we must convert A 7 = 
“TT = ~ A ^urm e an d ^7* = 1 Zm f rom natural units to SI units. We will use natural 

Lorentz-Heaviside units such that a = with eo = Mo = 1- To convert 

the magnetic field strength from natural Lorentz-Heaviside to SI units, we may 
consider for example the corresponding energy density in SI units: 


IB 2 

2 fJ,Q 


(A.l) 


We find that 1 Gauss (= 1CU 4 T) corresponds to 1.95 X 10 _2 eV 2 in natural 
Lorentz-Heaviside units. (See also footnote 24 of Raffelt and Stodolsky [36]). Using 
this conversion factor, we find: 


B 


A 7* = Sir = L53 >< 10 " 5 ( ^ 


2 M 


Bi 


10 14 GeV\ 

M ) 


kpc 


-l 


(A.2) 


We now turn to the conversion of A-, 


A 7 = 


— u: 


pi 


2 oj 


= 1.1 


Using a = |^ — ^ we obtain: 
n e \ ( keV x 


10" 2 cm -3 


u 


kpc 


-i 


(A.3) 


Comparing equations A.2| and |A.3 to the corresponding expressions in F2013 
(in between equations 9 and 10), we find that while our expressions for A 7 agree, 
F2013’s A 7 j is a factor of \[Ax\ too high, resulting in conversion probabilities that 
are a factor 4-7r too high. 


Appendix A.2. Treatment of photoelectric absorption 

Photoelectric absorption by the warm neutral medium is highly significant for 
photon energies in the ROSAT 1/4 keV band, but much less significant for the 3/4 
keV band. We include this effect using the standard density matrix formalism (see 
for example m) described in equations [l5] and [l8j As shown in figures [3] and [4j 


this leads to a substantial suppression of the ALP-photon conversion probability 
at c j = 200eV compared to that at u = 500eV (around 2 orders of magnitude). 

F2013 also considers photoelectric absorption, but uses a different propagation 
equation for the density matrix (see equations 10 and 11 in F2013). While he does 
not give his simulated conversion probabilities for energies in the 1/4 keV band, 
we can make some inferences from other information given. In section 2.2, F2013 
finds that similar conversion probabilities are needed to explain the 1/4 keV band 
and 3/4 keV band signals. However, his figure 5 shows similar M values required to 
explain the 1/4 keV and 3/4 keV bands. From this we conclude that his calculated 
conversion probabilities for the 1/4 keV and 3/4 keV energies are similar, and thus 
the effects of photoelectric absorption in the 1 /4 keV band have not been properly 
taken into account. 
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Appendix A.3. The ALP spectrum 
F2013 considers two ALP spectra: 


d$357 

dE 


~ 64.5 



E 

eV eW 


E 


v 

357 eVy 


cm 2 s 1 sr 1 eV 1 


(A.4) 


(his equation 4), and 


d &950 

dE 


~ 5.58 


AiVeff E 
0.5 eV 


exp 


E 

950 eV 


2 


cm 2 s 1 sr 1 eV 1 


(A.5) 


(his equation 5) 

In his equation 2, F2013 (correctly) states that these spectra should be nor¬ 
malised such that: 


AAeffPcMB = 5.9 x 10 2 AlV e ff eVcnr 3 (A.6) 

However, in fact both his spectra are normalised such that 

4 f°° 

Pcab = 2-7T- / ——EdE = 5.9 X 10~ 2 eVcm~ 3 (A.7) 

c J 0 dE 

i.e. he has presumably failed to include the factor of AN e g = 0.5. 

Furthermore, F2013 states that the energy of ALP spectrum in the 0.64-1 
keV band is 9.2 x 10 7 eV s _1 cm _2 sr _1 (the reported units are actually eV^ 1 , but I 
assume this is a typo). However, using F2013’s own spectrum, the energy is: 

[ EdE = 2.2 x 10 7 eVs _1 cm _2 sr _1 (A.8) 

40.64 keV dE 

Similarly, F2013 states that the energy of ALP spectrum d in the 100-284 
eV band is 9.3 x 10' eVs _1 cm 2 sr 1 . However, integrating the spectrum he gives 
in equation 4 we find 


PCAB = 


4 

4 \ 3 

IT 


pO.284 keV 


/ 0.1 keV 


d 4 > 357 

dE 


EdE = 2.3 x 10 7 eV s 1 cm 2 sr 1 


(A.9) 


These two mistakes combined mean that F2013 has overestimated the ALP flux 
itself by a factor of almost 10. We also note that F2013 uses an ALP spectrum with 
a higher average energy when attempting to reproduce the 3/4 keV band flux. Since 
the publication of F2013, CAB spectra with average ALP energies Eq.ab ^ 250 eV 
were excluded by overproduction of X-rays in galaxy clusters f8]. We therefore do 
not consider such models here. 
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Appendix A.4- Morphology 


F2013’s prediction of an observable X-ray signal from CAB to photon conversion 
in the Milky Way magnetic field is due to the compounding errors detailed above. 
We also note here that, as discussed in section 4.2, the morphologies of the Milky 
Way magnetic field and the unresolved X-ray background make a significant CAB 
contribution to the unresolved X-ray background impossible. This issue is not 
considered in F2013. 
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